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ABSTRACT 

Growing Supermassive Black Holes (SMBH) are believed to influence their parent 
galaxies in a negative way, terminating their growth by ejecting gas out before it 
could turn into stars. Here we present some of the most sophisticated SMBH feedback 
simulations to date showing that quasar's effects on galaxies are not always negative. 
We find that when the ambient shocked gas cools rapidly, the shocked gas is com- 
pressed into thin cold dense shells, filaments and clumps. Driving these high density 
features out is much more difficult than analytical models predict since dense fila- 
ments are resilient to the feedback. However, in this regime quasars have another way 
of affecting the host - by triggering a massive star formation burst in the cold gas by 
over-pressurising it. Under these conditions SMBHs actually accelerate star formation 
in the host, having a positive rather than negative effect on their host galaxies. The re- 
lationship between SMBH and galaxies is thus even more complex and symbiotic than 
currently believed. We also suggest that the instabilities found here may encourage 
the chaotic AGN feeding mode. 



1 INTRODUCTION 



Several properties of galaxies correl ate with the mass of 
SMBHs resid i ng in their centres (TTremaine et al.l 120021 ; 
lHaring fe Rixl 120041 ; ICattaneo etHf l2Q09h . suggesting a 
causal link. These relations are currently a ttributed to 
the apprecia ble destructive p owers of quasars (|Silk fe Reed 
Il998l : IZubovas fe Kin3 l2012h . At high SMBH accretion 
rates, quasars launch wide angle outflows with velo c- 
ity of v ou t ~ 0.1c as observed ([Pounds et al.l l2003bl lal; 
iTombesi et "all [2010) in AGN X-ray absorption spec- 
tra, ca pable of ejecting all the gas fro m the host 

galaxy ( Silk fc Reed Il998l; IZubovas fe King! l2012h . Ana- 
lytical arguments) King] l20Q3h set the momentum-driven 
critical SMBH mass limit, M a , that close l y matches the 
observed correlations jTremaine et al.l 120021 ; lHaring &; Rixl 
120041 : ICattaneo et al-lkood ). However how the outflow ex- 
pels the gas is not clear in detail: previous analyti- 
cal models include the essential outflow physics but as- 
sume spherical symmetry; numerical models are 3D but 
physically over-simplified due to computa t ional constraints 
(iDi Matteo et al.ll2005l : ISiiacki et aDl2007l : iBooth fe Schavd 
l2009h . 

As any fast outflow from a point source, the quasar out- 
flow leads to two shocks, one inner - the shocked SMBH out- 
flow itself - and the other outer, representing the shocked 
host galaxy gas. The two sho cks are separated by a con- 
tact discontinuity (see fig. 1 in lZubovas &; Kingll2012h . Pre- 
vious work has shown that the radiative cooling of the inner 
shock has a profound effect on the ability of the SMBH to 
drive the gas out of host galaxies. If the radiative time of 



the inner shock is short, only the momentum from the pri- 
mary quasar outflow is efficient in affecting t he host gas , 
with the mechanical energy lost to radiation ([Kind 2003). 
When the cooling time is long, the mechanical energy of 
the outflow thermalises (King 2005) and increases the hot 
gas pressure inside the quasar-driven bubble, leading to an 
even more powerful an d rapid v ~ 1000 km s _1 outflow 
([Zubovas fc Kindl2012h . 

Here we present numerical simulations that combine 
the strengths of the previous analytical (detailed feedback 
physics) and numerical (3D geometry) methods. Our simu- 
lations show that the way in which quasar outflows affect the 
gas in the host galaxies also strongly depends on whether the 
outer shock rapidly cools radiatively or not. In the former 
case, the shocked outer layer is very dense and is unsta- 
ble to the "thin shell" inst abilities previously known from 
supe rnova remnant studies ([Vishniad Il983l ; iMac Low et al.l 
1989). These instabilities lead to the shell breaking up into 
thin dense filaments that are resilient to the outward push 
from the quasar outflow but are susceptible to a triggered 
(or at least accelerated) star formation. We propose that the 
latter process may constitute a positive quasar feedback ef- 
fect on their host galaxies. Below we detail where and when 
this effect may be important. 



2 MODEL AND NUMERICAL METHODS 

Due to computational limitations, sub-grid models are a 
necessity to describe SMBH feedback in simulations that 
span cosmologically interesting volumes. Depositing energy 
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into the gas nearest to SMBH (e.g., SPH particle neigh- 
bours of the SMBH) is the most widely used method 
(iDi Matteo et al.ll2005l : ISiiacki et aDl2007l : booth fe Schavd 
2009) for such simulations. In this paper we study AGN 
feedback on scales of a single galactic bulge, and therefore 
we can afford (a) a much higher numerical resolution and 
(b) a more detailed physica l mod e l for AGN fe e dback. The 
outflow physics follows t he IKind d2003L l2005h: IKing et al.1 
(|201ll ) model (cf. §2.1 in IZubovas fe NavakshiiJl20ia ) while 
the numerical method follows IZubovas &; Navakshinl (|2012l ) 
exactly (see their §3.1), except for different initial condi- 
tions. In brief, hydrodynamics of gas in the host galaxy an d 
gravity solvers are based on code Gadget (|Sprin gel 2005). 
In contrast to the nearest neighbour method for AGN feed- 
back, we model the propagation of fast (v ou t = 0.1c) out- 
flow from th e SMBH explicit l y via "virtual particles" in- 
troduced by iNavakshin et al.l (|2009h . These particles are 
emitted isotropically by the SMBH, and propagate radi- 
ally outward with no self-interaction allowed. SPH density 
field is continuously calculated at the particles' instanta- 
neous locations, and the interaction with an SPH particle 
occurs only when the latter contains the virtual particle 
within its SPH smoothing kerne l. The method has been 
shown (INavakshin fe Powerl I20I0I ) to reproduce analytical 
results of IKind (|2003l 12005), and has been co mpared with 
the n e arest neighbours feedback m ethod by I Power et al.l 
(|201ll ). IZubovas fe Navakshinl (2012) further included opti- 
cally thin photo-ionisation and inverse Comp ton heating of 
the g alaxy's gas by the quasar radiation field ([Sazonov et al.l 
2005); this heating m ay be an important fee dback effect in 
gas-poor epochs (e.g., ICiotti &; Ostri ker 2007). 

Additionally, not only the momentum but also a frac- 
tion #e = exp[— Ri c /R] of the virtual particles' mechani- 
cal energy is passed to the ambient SPH particles, where 
Ric ~ 0.5 kpc (for the parameters chosen below) is the in- 
ner shock coolin g radius of the inner AGN outflow shock 
(|Kindl2003l . l2005l ). This radius is not to be confused with the 
outer shock cooling radius, R oc , introduced later in the Dis- 
cussion section of the paper. Note that g^(R <C Ric) — > 0, 
while gv(R ^> Ric) — ► 1- This energy deposition prescrip- 
tion mimics the transition from the radiatively efficient in- 
ner shock regime (small radii), when only the momentum 
of the SMBH outflow is passed to the gas, and the radia- 
tively inefficient regime (large radii), when both energy and 
momentum are availa ble to drive the host's gas ou t (see 
IKing et al.ll201lh . iFaucher-Giguere &; Quataertl (2012) have 
recently noted that if ions decouple thermally from electrons 
in the inner shock then the radiative cooling becomes far less 
efficient than obtained with a one-temperature model for the 
shock. The inner cooling radius Ri c in this case is negligibly 
small. We note that this (if correct) would introduce quanti- 
tative but not qualitative changes to our main results below, 
as we experimented with shells initially much smaller and 
also much larger than R oc (to be reported elsewhere). 

To demonstrate our points as clearly as possible, the 
host galaxy is modelled here by a fixed singular isothermal 
sphere potential softened at the centre where the SMBH is 
located. The ambient gas is initialised as a sphere with inner 
radius Ri n and outer radius R ou t = 10i?i n , at rest, and with 
the density profile given by 



p s (R) = f s po (R) = 



2ttGR 2 



(1) 



with p g (R) = for R < Ri n and R > i? ou t; here a — 
141 km s _1 is the velocity dispersion of the potential, f g — 
0.16/ is the gas ratio of the gas density to the underlying 
density of stars and dark matter; / > is a free parameter. 
Initial gas temperature is set such that the sound speed is 
equal to a. At time t = 0, the quasar outflow is turned on. 
The critical M a mass above which the SM BH outflow should 
expel the gas for this potential (Navaks hin &; Powerl [2010) 
is M a « 1.5 x 1O 8 /M . 

We fix the quasar momentum outflow rate and the 
SMBH mass in order to avoid model-dependent compli- 
cations near SMBH. To avoid increasingly short time- 
steps at small rad ii, we use the accretion radius approach 
(Po wer et al.ll201lh and remove SPH particles from the sim- 
ulation that are closer than 30 pc from the SMBH. The 
densest and coldest self-gravitationally bound regions would 
continue to contract to arbitrarily large densities, a process 
which physically terminates in star formation. When gas 
density exceeds 200 times the potential's tidal density po, 
and the Jeans mass in the clump is smaller than the SPH ker- 
nel (40 particle neighbours) mass, SPH particles are turned 
into collisionless star particles. Our treatment of star forma- 
tion neglects feedback effects from massive stars; this does 
not alter our main conclusions. 



3 RESULTS FOR SPHERICAL SHELLS 

We first present two simulations that exemplify the effects of 
the quasar outflow on the host galaxy's gas in two extremes, 
e.g., (i) when the shocked gas can cool rapidly and (ii) when 
radiative cooling is inefficient (low gas density). 



3.1 Thin shell instabilities 

Figure [1] shows the density (left column) and temperature 
(right column) of the ambient gas shocked by the outflow 
in simulation Fl in the gas-rich epoch (/ = 1, top panels), 
and in simulation F0.03, corresponding to a gas-poor epoch 
regime (/ = 0.03, bottom panels). For both simulations, the 
SMBH mass is M D h — 10 8 M and the outer shell radius 
is -Rout = 2 kpc. The host galaxy is rapidly emptied of gas 
in F0.03; the shock front velocity is about 1500 km s _1 , so 
that the shock expanded to R ~ 1 kpc by the time of the 
snapshot, t — 0.57 Myr. In contrast to that, the same SMBH 
outflow in the gas-rich simulation Fl stalls: the radius of 
the shock front is R ~ 0.4 kpc at time t = 3.8 Myr. This 
immense difference is due to two factors, both significant 
here: (i) the weight of the gas larger by a factor of « 30 
in Fl than in F0.03, and (ii) the inner shock remains in 
the momentum driven regime in Fl, since the shell stalls 
at R < Roc, whereas the shock becomes non-radiative and 
thus energy driven in simulation F0.03. These results are 
consistent with t he analytical expectat ions for spherically 
symmetric shells ([Zubovas &; Kindl2012h . 

An unexpected result is that while the shocked shell is 
quasi-spherically symmetric in F0.03, it goes violently un- 
stable in Fl. In hindsight, the result is completely analo- 
gous to the stability of supernova remnant shells. Cold dense 
swept-up shells are strongly unstable to Ray leigh- T aylor and 
IVishniad (| 19831 ) instabilities (|Mac Low et al.l ll989), whereas 
non-radiative hot and geometrically thick shocks are stable 
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(compare fig. 3 in iMac Low et al.l (p.989) with the rest of 
their figures). Physically, radiative cooling is key for devel- 
opment of the instabilities because any initial density in- 
homogeneities are amplified when denser regions cool and 
get compressed by the high pressure surroundings; without 
cooling these density inhomogeneities tend to diffuse away. 

The cold dense filaments are crucial to understanding of 
SMBH feedback on host galaxies and also the SMBH feed- 
ing, as we shall see. 

3.2 Accelerated star formation 

The cold dense filaments in simulation Fl are the birth place 
for newly born stars. The right panel of Figure [2] presents 
the gas temperature projected along rays as seen from the 
SMBH for simulation Fl. The shell is a web of cold dense 
filaments with the strongest density peaks occurring at fil- 
ament intersections. The left panel of the figure shows the 
density and velocity field at the time when dense filaments 
are starting to fall inward. Stars (red or cyan dots) are born 
inside the densest and coldest filaments. 

Physically, quasar shocks are able to trigger star forma- 
tion, or accelerate it, by compressing the cold gas in the 
hosts because the pressure produced by the quasar out- 
flow is far higher than the maximum gas pressure inside 
the pre-quasar host, Pmax ~ fgpocr 2 (if the gas were hotter, 
c 2 > a 2 , a thermal pressure-driven wind would result). The 
pressure in the shocked ambient gas and inside the hot bub- 
ble inflated by the quasar is P s h ~ fgpov 2 , where v e ^> o 
is the energy-dri ven shell expansion velocity calculated in 
iKing et alj feoilh . In the energy-driven regime, the quasar's 
bubble pressure is higher than P ma x by the ratio 

Psh 
P d 

To exemplify the quasar's ability to trigger star formation 
further, consider the effects of the quasar outflow on a star- 
burst gas disc. Such discs may have star formation rates 
as high as ~ 1 s yr -1 while self - limiting their star 
formation rates (Thomps on et al.1 [2005) by star formation 
feedback. Their midplane pressure is Pd ~ poo~ 2 fd, where 
fd < fg is the disc mass fraction at a given R. This is smaller 
than Pmax , so we see that even the starburst-supported discs 
are compressed to higher densities by quasar outflows. Since 
star formation rates ar e proportional to Pd in these discs 
([Thompson et a 

D l2005h , disc compression due to quasar 
shock results in an even faster starburst in the disc. 



3.3 Chaos and order in SMBH accretion 

SMBH accretion is still far from understood. One of 
the significant challenges is the "star for mation catastro- 
phe" of self- grav i tating parsec scale discs (|Paczvhskilll978l : 
Nayakshin et al. 2007) in which gas is consumed more 
rapidly in star formation than it is fed t o the SMBH. On e 
proposed solution is disc self-regulation ([Goodman! 12003 ), 
and the other is "chaotic accretion", in which angular mo- 
mentum of g as is cancelled in sho cks rather than trans- 
ferred away ( King fe Prinde] l200d ). Chaotic SMBH feed- 
ing is suggested to enco urage a more rapid S MBH growth 
in the early Universe ([King &; Pring le 2006) by limitin 
SMBH spin and may also explain ( Navakshin et al1l201 



24a-o /3 /- 2/3 



> 1 . 



(2) 



why SM BH do not correlate we ll with properties of galac- 
tic discs ( Ko rmendv et al.l l201lh . Such flows could develop 
as a result of messy ini tial conditions or s upernova- driven 
turbulence in the bulge (Hobb s et al. I l201lh. 

Thin shell instabilities (Vishniac 1983), re-discovered 
here in the SMBH growth context, are thus an additional 
way of developing such multi-phase chaotic inflows start- 
ing from non-turbulent spherically symmetric initial condi- 
tions. The resulting gas flows are indeed highly variable as 
infalling fingers (see figure 2) shadow the ambient gas in an 
unpredictable way. On the other hand, if the ambient gas 
possesses symmetry before SMBH feedback is turned on, 
then there may be preferred inflow directions and an orderly 
flow/ou tflow develops in this case: fo r example, see simula- 
tions bv lZubovas &; Navakshin] ([2012) for the Fermi Bubbles 
in the centre of the Milky Way. Further simulations with a 
range of realistic initial conditions are needed to investigate 
the importance of chaotic inflows for SMBH feeding. 



4 IMPLICATION FOR NON-SPHERICAL 
GEOMETRIES 

Having investigated an initially spherically symmetric shell 
in 33 we now wish to ask what happens in a more realistic 
case. Based on the earlier results, we expect that deviations 
from spherical symmetry in initial conditions will be ampli- 
fied by the thin shell instabilities when the ambient shocked 
gas is able to cool rapidly. As before, the flow may break into 
dense inflowing and lower density outflowing regions or di- 
rections. If angular momentum is important then the denser 
flows may settle into a circularised disc. 

We wish to emphasise the fact that the cooling rate 
of the shocked ambient gas is key to this behaviour. Fig- 
ure [3] presents two contrasting simulations with same ini- 
tial conditions but with ambient gas radiative cooling either 
on (left panel) or off (middle panel). The initial conditions 
are obtained from the spherical shell used in Fl by a linear 
transformation z —> z/2 which obviously makes an elliptical 
rather than a spherical initial shell. The shell also slowly 
rotates, v ro t = 0.3v c irc, where v c irc = 2 1 ^ 2 a is the circular 
velocity for the potential, with the angular momentum vec- 
tor pointing in the positive z direction. Finally, the SMBH 
mass is slightly increased to Mbh = 1.5 x 1O 8 M0, the M a 
mass for this potential. 

The left panel of Figure [3] shows that the outflow is 
channelled along the symmetry axis. However, as the ambi- 
ent gas cools, most of it collapses onto the symmetry plane, 
becoming resilient to feedback and instead mainly turning 
into stars. This is in a stark contrast to a simulation shown 
in the middle panel of the same figure, where the radiative 
cooling in the ambient gas is artificially turned off. All of 
the gas is expelled in this case because the shocked ambient 
gas remains hot and always presents a large enough surface 
area for the SMBH feedback to work on. 

To quantify this cooling-regulated behaviour, we intro- 
duced a free parameter A coo \ ^ 1, which is used to artificially 
reduce the radiative cooling rate function A — >> A/A coo i in a 
set of numerical experiments which share same initial con- 
ditions as explained above. The left panel of Figure [3] cor- 
responds to A coo \ = 1, e.g., the proper unmodified cooling 
rate, whereas the middle panel of the figure corresponds to 
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A CO o\ = oo (e.g., cooling in the ambient gas is turned off). 
The right panel of figure [3] shows the fractions of gas ex- 
pelled and turned into stars as functions of A coo \. We see 
that quasar outflows mainly turn gas into stars when the 
shocked gas cools rapidly but is able to expel most or even 
all of the gas if the shocked gas is unable to cool (large 



5 DISCUSSION 

We presensted simulations of the interaction between the 
fast outflow emanating from a quasar and the ambient gas 
in the host galaxy utilising the he virtual particle method 
of iNavakshin et al.l (|2009l). We believe t he method is an 
improvement (see also I Power et al.l 1201 lh on the sub-grid 
nearest-particle methods since the SMBH outflow is explic- 
itly modelled to track where and how it impacts the am- 
bient gas, and whether the outflow's momentum, energy, 
or both are passed to the ambient gas. In particular, con- 
sider the left panel of figure [2] The cold filaments subtend 
a small solid angle as seen from the quasar location and 
should thus encounter a small fraction of the quasar outflow 
only, with most quasar feedback impacting on the low den- 
sity gas instead. This is correctly captured by our virtual 
particle method. In contrast, the nearest neighbour meth- 
ods may spread all of the quasar feedback energy to the cold 
filaments because they tend to be closer to the SMBH. 

Our key new result is the finding that the response of 
the ambient gas to the fast outflow of gas from the quasar 
depends critically on whether the ambient gas is able to 
cool rapidly or not. Generally, in the latter case the ana- 
lytical criteria, such as the outflow momentum (thrust) bal- 
anced against the weight of the gas in the bulge are suffi- 
cient to establ ish whether the s hell is driven out ward, stalls 
or collapses (|Kindl200l 120051 : king et alJ l201lh . However, 
when the outer shock is radiative, even initially spherically 
symmetric shells develop a multi-phase structure. Denser 
regions cannot be driven outward easily. Quasar outflows 
affect these high density regions most effectively not by ex- 
pelling them outward but by compressing them and trigger- 
ing (or at least accelerating) star formation there. 

Let us now quantify the conditions when the ambient 
shocked gas is able to cool rapidly. Bremsstrahlung is the 
most significant cooling function for the problem at hand, 
with A(T) = A T 1/2 , A = 3.8 x 10" 27 erg cm 3 s _1 . The 
shocked gas cooling time, t coo \ — 3kT/nA(T), where n and 
T is the shocked gas electron density and temperature, re- 
spectively. In the strong shock limit, n — 4p g (R) / '(firrip), 
where \i — 0.63 is the mean molecular weight and p g (R) 
is the pre-shock gas density given by equation [T] and T = 
(3/16)/j J m p Ve/kB- To determine whether the quasar shock 
is radiative or adiabatic, t coo i should be compared with the 
shocked shell flow time scale, R/R, where R = (4/3) v e is 
the outer shock velocity (| Zubovas &; Kind 1 20121 ): 



D 1 Ol ^»5/3 2/3 / M a \ 2/3 



Ac 



R/R 



0.83 3-8 X 10-" B/ , e 
Ao 
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(M hh \ 

\ M a ) 
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#k P c • (3) 



At a fixed fraction /, the "outer cooling radius", R oc , the ra- 
dius which separates the rapidly cooling outer shock regime 
from the slowly cooling one, can be defined by requiring that 
the ratio in equation [3] be equal to one. The result is: 



(4) 

3.8 x 10- 27 v ; 

In the early gas rich epochs (which we operationally de- 
fine as / > 1) SMBHs affect their hosts both positively via 
accelerating star bursts at "small" radii (at R R oc ), and 
negatively at larger radii where the gas is more likely to be 
blown out. Interestingly, the powerful starbursts in submil- 
limeter galaxies (e.g., galaxies forming stars at rates as high 
as ~ 10 s Mq yr 1 a t redshift z ~ 2) have been resolved 
bv lTacconi et all I (|2006h to have intrinsic sizes of a few kpc, 
comparable to R oc in the gas-rich epoch. 

In the later gas-poor epochs, when / < 1 everywhere, 
and Mbh ^ M a , it follows from equation 2]t hat R oc <C 1 kpc, 
so that shocked gas is unlikely to clump up into high density 
features over most of the observationally interesting scales 
in the host galaxy. In this regime quasars should be better 
able to sweep the hosts clear of gas if they are activated 
for long enough. This is the widely appreciated negative 
quasar/ AGN feedback mode, and it probably explains why 
the present day ellipticals ar e "red and dead" and w hy their 
SMBHs are not very bright (|Schawinski et alJl2009h . 



6 CONCLUSIONS 

Our results broaden the spectrum of possible effects that 
SMBHs have on their host galaxies, showing that it is not all 
about ejecting the gas out. We found that it is increasingly 
difficult to drive gas out of the hosts at small radii, R < 
R oc (equation [4} , in the early gas-rich epoch. Even in an 
initial spherical symmetry the thin and dense outer shell 
breaks into filaments and clumps that are hard to drive away 
but which can be readily over-pressurised by the hot quasar 
wind and which result in enhanced star formation rates. The 
latter process is a form of positive rather than negative AGN 
feedback on host gala xies (we note that sim ilar statements 
were already made by ISilk fc Normap||2009L based on their 
analytical arguments) . Non-spherical initial shell geometries 
only strengthen these conclusions as the multi-phase inflow- 
outflow separation becomes even more robust (see 2}. 

We also found that in later gas-poor epochs, when the 
shocked host's gas cannot cool rapidly (cf. Discussion), the 
inhomogeneities do not develop (cf. the bottom panel of fig. 
1) or do not grow for non-spherical initial conditions (cf. 
figure 3, the middle panel). The SMBH is then an efficient 
negative feedback engine as has been appreciated for a long 
time now. 
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Figure 1. Gas outflows for a gas-rich simulation Fl (top panels) at time t = 3.8 Myr (top) and same for a gas-poor simulation F0.03 
at t = 0.57 Myr (bottom panels). The left panels show projected gas density and the right ones show gas temperature. To expose the 
3D nature of the flows, the left sides of each panel show geometrically thin "wedge-slice" projections of the gas flow, where only regions 
y 2 /(x 2 + y 2 + z 2 ) < 1/16 are shown; the right-hand sides of the panels sho w full cube projections. The shell in the gas-rich simulation 
cools rapidly and becomes unstable to thin shell instabilities ( Vishniac 1983). The shell in the gas-poor simulation cools very slowly, and 
is thus hot and radially extended. This shell is stable and remains spherically symmetric as the gas is driven away. 
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x [ kpc ] <p I radian ] 

Figure 2. Left panel: "Wedge projection" (see fig. 1) of gas density and velocity vectors for simulation Fl at time t = 5.2 Myr. The 
high density regions of the shell receive less quasar feedback per unit mass and thus fall inward. Red symbols show newly born stars 
that form within the dense filaments. Most of the stars form in clusters which are disrupted when passing through the central regions 
of the galaxy. Right panel: Temperature of the gas as seen from the SMBH, averaged along rays, for a shell 0.2 < r < 0.4 kpc, at time 
t = 4.3 Myr, and as a function of the two spherical angles in spherical coordinates, defined in the usual way: cos# = z/ sj x 2 + y 2 + z 2 , 

and sin</> = y/ sj x 2 + y 2 . Cyan symbols show locations of the stars. These form exclusively inside the density peaks (corresponding to 
minima in temperature), which occur most frequently at the intersections of the filaments. The colour scales for the panels are same as 
they are in figure 1 for simulation Fl. 




Figure 3. Quasar feedback on an initially elliptically-shaped and slowly rotating shell in the cases with the proper radiative cooling 
on for the shocked gas (left panel) and for the case when radiative cooling is switched off (middle panel). In both cases gas is expelled 
along the vertical direction. However, with cooling on, the gas settles into a high density disc that is resilient to quasar feedback. Very 
little gas is ejected from the galaxy in this case; most is turned into stars or accreted through the inner boundary condition. In contrast, 
with cooling off, all of the ambient gas is expelled to infinity because no high density features form; hot and diffuse gas is much more 
susceptible to quasar feedback. The right panel shows how the fraction of the shell expelled to infinity and turned to stars at the end 
of the simulations (t = 9 Myr) change as radiative cooling is progressively suppressed by factor A coo i (the left panel corresponds to 
^-cool — 1> an d the middle panel to A cool = oo). More and more gas is ejected rather than turned into stars as cooling is suppressed. 
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